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COMPUTATION OF GENERAL PLANETARY PERTURBATIONS 
FOR RESONANCE CASES 


INTRODUCTION 

Many investigators have studied resonant motions in celestial mechanics. 
These problems are of theoretical interest from the point of view of stability 
and for understanding the characteristics of the motion. The practical impor- 
tance of these problems comes from the existence of so many cases of resonant 
motions in the solar system. Hagihara's work (Reference 3) contains an ex- 
cellent summary of studies of these problems. The work of Schubart (Reference 
5) has added significantly to the understanding of motions of resonant asteroids. 
By a numerical averaging process on a computer he was able to isolate the 
secular and critical terms of the disturbing function and then investigate the 
stability and characteristics of the motion for a wide range of values of the 
parameters. 

The cases of resonance require special techniques when general perturba- 
tions are to be computed. For this purpose it is necessary to represent the 
motion in a form compatible with the usual trigonometric series so that accurate 
results may be obtained including higher order effects from all the disturbing 
planets. Such a form is obtained in the present study by using an argument in the 
trigonometric terms such that the resonance of the mean motions is taken into 
account. In this way the critical terms are removed and the effects appear in 
the secular terms. The resulting series have the same form as in ordinary 
nonresonance cases. The mean anomalies of the two bodies have been replaced 
by a single related "anomaly" in the arguments of the trigonometric terms. 


THE REFERENCE ORBITS 

The perturbations of the coordinates are computed to give the deviation of 
the motion of a planet from the motion in a fixed Keplerian ellipse. There is 
some latitude in the choice of the elements of the reference ellipse, but to keep 
the perturbations small, the reference motion should be a close approximation 
to the actual motion. When the mean motions, n and n' , of two planets are in 
the ratio of two integers, 


n p 

— = — , p, qare integers, 
n q 
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the motions are said to be commensurable and this is a case of resonance which 
requires some special consideration in the computation of general perturbations. 

An extensive expansion of the disturbing forces from the mutual attractions 
of two planets is given in Reference 2. In the usual nonresonance case the per- 
turbations are obtained as trigonometric series in terms of the two mean anom- 
alies, g and g' , with various powers of the time as factors. This form results 
from the fact that the position vectors ? 0 and in the two reference orbits are 
periodic functions of the mean anomalies g and g' respectively. In the expan- 
sions of the components of the disturbing force, the terms in the series have the 
form 


where i and j are integers. When these terms are integrated with respect to 
the time, the coeffecients are multiplied by the factor 


n 

in + jn' 


In cases of resonance where 


qn = pn' 


as before the terms with the argument 


qg - pg' 

or any multiple of this become critical since 


qn - pn' = 0. 


There are several possibilities for treating such cases. The purpose here is 
to give a method wherein the form of the series is similar to that of the usual 
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nonresonance case. The mean anomalies, g and g' , are linear functions of the 
time, t : 

g = nt + g Q 
g' = n't + g Q '. 

Putting 

n =— = — and g* = n*t 
P q 


there follows 


g = pn*t + g Q = pg* + g 0 
g' = qn* t + g ; = qg* + g ' 

Now if g*is increased by 2 w theng and g' are increased by 2pw and 2qw respec- 
tively. Since r 0 and r' Q are periodic functions of g and g' respectively, both 
r (j and ?' are periodic functions of g*. The former development in terms of 
g and g' can be replaced, in the resonance case, by a development in terms of 
the single angle g*, and in this form there are no critical terms. In the results 
the affects of the critical terms will appear in secular form. The development 
given on pages 5 through 11 of Reference 2 can be used in the computation ex- 
cept that the double harmonic analysis is now replaced by simple harmonic 
analysis in the variable g*. 

In the application of this method, the mean motions of the two planets in the 
reference orbits must be exactly commensurable. The deviations from the 
reference orbits are properly taken into account by adjusting the constants of 
integration in the series. 


FORMULAS FOR THE PERTURBATIONS 

The form presented in the previous section for representing the motion in 
cases of resonance using trigonometric series can be applied using nearly any 
decomposition of the perturbations. Musen’s (Reference 4) method of pertur- 
bations of the rectangular coordinates is used in the present case. 


3 



The position vector, r , of a planet is given by 


7 = (1 + a) 7 Q + /3w + yaR 


where a, /3 ,y are the components of the perturbations, r Q is the position vector 
in the fixed reference ellipse, a is the semi-major axis of the reference ellipse, 
R is the unit vector in the direction of the angular momentum of the motion in 
the reference ellipse and 



w — , 

n dt 

n being the mean motion in the reference ellipse and t is the time. 
The perturbations are computed from the formulas 


a = J (MjD + M 2 E) d (nt) 
P = f (M 4 - 2a) d (nt) 


7 - 


J M 3 Dd(nt) 


to which the contributions from the constants of integration are to be added. 
The integrands are obtained from the formulas: 


Mj = Sj ‘ (a 2 F) 
M 2 = S 2 • (a 2 F) 

M 3 = ~ R ' (a 2 F) 
M 4 = jM 2 d (nt) 
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sin e 


S. = P cos e + 0 

1 /TT? 


s 


2 


- P sin e + Q 


/f~- e 2 


and 



[sin (r] - e) 


e sin 7 } + e sin e] 


E = 2 — [1 - cos (17 - e) ] 
r o 

where e is the eccentricity of the reference ellipse, e is the eccentric anomaly 
in the reference ellipse, 77 is the eccentric anomaly held constant during the 
integration and replaced by e afterward, and P ,Q are the standard unit vectors 
associated with R . The factor a 2 F is the disturbing force which can be ex- 
pressed as 


a 2 F = a 2 f j + a 2 f j t + a 2 f 2 t 2 + 


The various f ’s which are factors of the powers of the time are given on pages 
6-10 of Reference 2. 

The f 's depend on the reference position vectors of r Q and r' Q which are 
periodic functions of g*. Similarly S x , § 2 , and r Q are periodic functions of g*. 
Thus , M 2 , and M 3 can be written in the form 


M j = lUj i + lTlj jt + tllj 2 t 2 + • * * 

where 

ftl . = S. ‘a 2 f. 

1 . j 1 j 
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m 2i . = s 2 • a*f. 


m . - — - R • a 2 f . 
3 .J a J 


for j = I, 1,2, 3 Therefore, by simple harmonic analysis each Hlj } 

can be obtained as a trigonometric series in g*. 

The expressions for D and E also can be represented as trigonometric 
series in g* and 77. These series along with those for the Mj's form the in- 
tegrands in the expressions for a, fi, and y. Terms under the integral have 
the general form 


[C cos (ax + b) + S sin (ax + b)] x p . 


Upon integration this takes the form 


[ cos (ax + b) + S p sin (ax + b)J x p + £ C p _ ^ cos (ax + b) + S p _ j sin (ax + b)J x p 

+ [ C p _ 2 cos (ax + b) + S p _ 2 sin (ax + b)J x p_2 + • • • 

+ [C 0 cos (ax + b) + S 0 sin (ax + b)J 
where 


S 


a 


C 

* 

a 


C P-i=aV S p-i-T C p 


q - R \ s 

V2 - n Vi ■ 


p- 1 


p- 2 


“p- 1 ’ 


Q - R S 

p - 3 a V2’ 


P-2 


P-3 


~P“2 ’ 
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and finally 


C 


o 




1 • 


The special case where a = 0 must be treated separately. Here the term 
under the integral has the general form kt p where k is a constant. Integration 
yields terms of the form(k/p + l)t p+1 which introduces factors with one higher 
power of the time and since /? involved double integration factors are generated 
for two higher powers of the time. In this manner, secular terms and terms 
with higher powers of the time enter the expansions. 


The final expressions for a , /3, and y have the general form 


£ t p £[ C P .i cos (ig*) + S pii sin(ig*>] 

p = 0 i L 

where s is the highest power of the time carried in the computation and g* = n*t. 


APPLICATION 

The case of the mutual attractions of Pluto and Neptune demonstrates the 
capabilities of the method. The values used for the reference orbit elements 
are given in Table I. These are approximate mean elements with the mean 
motion and semi-major axis of Pluto adjusted to satisfy the resonance criterion, 
3n g = 2 n g where n g is the mean motion of Pluto and n g is the mean motion of 
Neptune. 


Table I 

Reference Elements of the Planets Pluto and Neptune 


Planet 

Semi -Major 
Axis in a.u. 
a 

Eccentricity 

e 

Mean Motion 
in 

Degrees/Day 

n 

Longitude 

of 

Ascending Node 

Q 

Argument 

of 

Perihelion 

u 

Inclination 

i 

Mean Anomaly 
at Epoch 
(JD 2415200.5) 
g o 

Neptune 

30.005162 

.008956 

?0059818842 

131?2332 

275?9147 

1^7745 

39?1226 

Pluto 

39.382760 

.248895 

! 0039879228 

109?6750 

1 13?9034 

17^1434 

230?0159 
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The coefficients of T k cos(ig*) and T k sin (ig*) in the series for the 
perturbations are given in Tables II and HI. For comparison purposes, the 
mutual attractions of Pluto and Neptune were computed inChebyshev series, as 
in Reference 1. The series in Table II and the Chebyshev series were evaluated 
at intervals of 800 days to give residuals in a , /3 and 7 for two centuries on either 
side of the epoch. Figure 1 shows the trend of the residuals over the interval. 
The unit of time for the secular terms in the trigonometric series is one century 
(36525 days) and as expected the residuals increase in size beyond one unit from 
the epoch. The largest deviation of the four centuries is .7 x 10 ' 8 or less than 
002 of arc in (S. 


Table II 

Coefficients in the Series for the Perturbations of Pluto Due to Neptune 


0.1O 6 


i 

j 

k 

cos 

sin 

cos 

sin 

COS 

sin 

0 

0 

0 

7. 7108 

0. 

0. 

0. 

3.0499 

0. 

1 

0 

0 

116.4417 

76.7693 

237.3396 

-250.1 348 

10.2250 

-4.3397 

2 

0 

0 

20.4367 

-6.8767 

- 0. 0000 

O .0000 

- 0 . 0 0 p. n 

-P . POPP 

3 

0 

0 

13.8694 

-10.2346 

-5.1872 

-12.4812 

1.1377 

-6.1786 

4 

0 

0 

7 • 8 8 1 2 

-3.0165 

-7.3918 

-1 1.2340 

-2.2129 

2.6157 

5 

u 

u 

4 . 0 0 0 2 

0.0170 

-1.3823 

-6.214^ 

-0.6797 

0.6806 

6 

0 

0 

1.4731 

-0.3964 

-0 . 6360 

-2.3503 

-0.2082 

0.2638 

7 

0 

0 

0.3957 

-0.2325 

-0.4682 

-0.7887 

-0.1218 

0.3940 

8 

0 

0 

0.4234 

-0.2927 

-0.4400 

-0.4676 

-0.C971 

0.2566 

9 

0 

0 

0.2797 

-0.2631 

-0 . 3533 

-0.31 40 

-0.0672 

0. 1060 

10 

0 

0 

0.1484 

-0. 1004 

-0. 194 0 

-0. 1977 

-0.0233 

0.0645 

11 

0 

0 

0.0969 

- 0 .05 9 8 

-0. 1064 

-0.1157 

0.0002 

0.0557 

12 

(} 

0 

'‘.0501 

-0.0554 

-0.0735 

-0.0573 

-0.0027 

0.0398 

1 3 

0 

0 

n.0262 

-0.0368 

-0.0556 

-0.0305 

-0.0054 

0.0243 

14 

0 

0 

u.0191 

-0.024 1 

-0.0388 

-0.0209 

-0.OO12 

0.0141 

15 

0 

0 

0.0116 

-0.0180 

-0.0257 

-0.0134 

0 . .0 0 1 6 

0.0093 

16 

0 

0 

0.0062 

-0.0121 

-0.0167 

-0.0071 

0.0014 

0.0072 

17 

0 

0 

0,0041 

-0.0080 

-0.0116 

-0.0036 

0 • 0 0 0 7 

0.005 1 

18 

0 

0 

0.0024 

-0.0057 

-0.0083 

- 0 .0020 

O • n 0 0 6 

0.0031 

19 

0 

0 

4.001? 

-O.0039 

-0.0067 

-o.ooll 

0.0O07 

0.0019 

20 

0 

0 

o.onofi 

-^.00?7 

-0 . pp?8 

-n . OOOf) 

0 . n r 0 7 

P . P p 1 4 

0 

0 

1 

92.0432 

0. 

0. 

0. 

1 .3160 

0. 

1 

0 

1 

-0.0945 

-0.0592 

-0.868C 

0.7988 

-0.0O07 

0.0034 

2 

r 

1 

“14.0529 

-8 1.7793 

-163.6033 

28.0364 

10.0917 

3.2832 

3 

0 

1 

-0.0483 

-0.000b 

0.0381 

0.0756 

-0.0024 

-0.002 1 

4 

G 

1 

8. mi 

3.1016 

6. 104 1 

- c . 7 6 4 5 

-1.1030 

0.6854 

5 

r; 

1 

0.0046 

0.0016 

— 0 .0032 

0.0037 

0 . 0 0 ] ? 

-0.0016 

fa 

0 

1 

-1. /391 

0.6 39 7 

0 • h 2 7 4 

1 • 1 7 -> 2 

0*0345 

- G • ^ 3 7 6 

7 

0 

1 

—0 .004 0 

•>.0014 

0 .0014 

O.0011 

0.0002 

-0.0001 

8 

0 

1 

0.1385 

-0.3828 

-0. 1910 

-0.0664 

0.0349 

0.0387 

0 

0 

1 

0.0002 

-0 . 0 0 no 

0.0042 

0 . 0 n n 7 

0 . n 0 r, n 

-p . 0002 

10 

{) 

1 

0.0486 

0 • ( 1 0 4 6 

0 . 0 ? 3 9 

-0.0192 

-0.0124 

0.0003 

11 

( ; 

1 

-0.00O1 

r .00^5 

?.orr 1 

- .0002 

p . PPP 1 

-n .non 1 

12 

0 

1 

-0.0246 

—0.004 2 

- 4 . CO 1 3 

^ .0 08 

0.0020 

-0.0025 

13 

0 

1 

-0.0002 

- 0 .0 00'' 

o.o 0 0 1 

4 .000 1 

0 • p p 0 0 

- 0 . 0 0 n 0 

14 

0 

1 

0.0050 

-0 .0042 

-0.0012 

-0.0 014 

P.0002 

0*0008 

16 

0 

1 

- 0.0000 

n.OGin 

0 . 0 n 0 5 

0.0 0 0 '*■ 

- P . P P P 2 

-p. OOP 1 

18 

0 

1 

-0.00O4 

- 0 .onn 8 

- 0.0001 

n.onni 

p . 0 r 0 1 

- 0 . pppp 

20 


1 

• 0 0 0 1 

— n . r- 0 0 4 


_r # n'' A r 


p # ppp p 
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Table II (Continued) 

Coefficients in the Series for the Perturbations of Pluto Due to Neptune 


tt.10 6 0 .io 6 r.io 6 


i 

j 

k 

COS 

sin 

cos 

sin 

cos 

sin 

0 

0 

2 

0.0191 

0 • 

-188.5243 

0. 

0.0010 

0. 

1 

0 

2 

. 

CD 

-0.1433 

. 

0 

1 

-0.4173 

0.0059 

-0.0025 

2 

0 

2 

0.0390 ■ 

0.0329 

0.0700 

-0.0767 

0 • 0 0 n 8 

-0.0025 

3 

0 

2 

0.0017 

0.0206 

-0.0003 

0.0093 

-0.0005 

-0.0018 

4 

0 

2 

-0.0144 

0.0004 

-0.0032 

0.008 3 

0.0012 

0.0005 

3 

0 

2 

-0.0003 

-0.0046 

-0.0021 

0.0032 

0*0006 

0.0003 

6 

0 

2 

0.0008 

-0. 0048 

-0.0014 

0.0013 

-0.0001 

0.0000 

7 

0 

2 

0.0020 

-0.0018 

-0.0008 

n . n no 7 

O.oooi 

-0.0000 

8 

o 

2 

0.0004 

0.0011 

-0.0004 

0 . 0 00 5 

0 * ono? 

-o.oooi 

9 

0 

2 

-0.0003 

o.OOOB 

-0.0001 

0 . OOn 3 

O.oooi 

-0 .onon 

10 

0 

2 

-0.0007 

-0.0002 

-0.0001 

0.0002 

0.0001 

0.0000 

11 

0 

2 

-0.0003 

-0.0003 

-O.0001 

n . 0 002 

o.oooo 

0.0000 

12 

o 

2 

0.0001 

-0.0002 

-0.0001 

0 .000 1 

0.0000 

-o.oooo 

13 

0 

2 

o.oooi 

-0.0000 

-0.0000 

o.oooi 

0.0000 

-0.0000 

14 

0 

2 

-0.0000 

0.0001 

-0.0000 

0.000 1 

o .oooo 

-0 . onno 

15 

0 

2 

-0.0001 

0.0000 

-0.00,00 

o.oooo 

O.OOOO 

-0.0000 

0 

0 

3 

-0.0122 

0. 

-0.0345 

0. 

-0.0005 

0 . 

1 

0 

3 

o.oooi 

0.0002 

0.0008 

-0.0007 

0.0000 

-o.oooo 

2 

0 

3 

-0.0298 

0.0179 

0.0172 

0.0055 

-0.0014 

0.0029 

3 

n 

3 

-0.0000 

c: 

c 

c 

. 

c 

1 

0 . 0000 

-O.OOOI 

n . noon 

o.oooo 

4 

0 

3 

0.0041 

-0.0123 

-0.0001 

0.0011 

-0.0003 

-0.0006 

6 

0 

3 

0.0024 

0.0044 

-0.0002 

-0.0001 

0*0002 

0.0000 

8 

0 

3 

-0.00] 8 

-0.0004 

0.0000 

-0.0000 

-0 .oooo 

0.0000 

10 

0 

3 

0.0005 

-0.0004 

o.oooo 

o.oooo 

n .0000 

-o.oooo 

12 

f) 

3 

-0.0000 

0.0002 

-0.0000 

-o.oooo 

O.OnOO 

0.0000 

14 

o 

3 

-O.0OO] 

-o.oooi 

0.0000 

-n „ noon 

-0. on. 0 0 

n.0000 

0 

G 

4 

-0.0196 

0. 

0.0227 

0. 

-0.0000 

0. 

1 

0 

4 

-o.oooi 

-0.0001 

0.0001 

0.0OO2 

-0.0000 

0.0000 

2 

n 

/. 

0-0092 

-0.0110 

-0.0000 

0.0000 

-0.0000 

0.0000 

3 

0 

4 

0.0001 

o.oooo 

0 . 0000 

-0 .0000 

0.0000 

0.0000 

4 

0 

4 

0.0009 

0.0051 

0.0000 

-0.0000- 

-0.0000 

0 • 0 o o o 

6 

0 

4 

-0.0016 

-0.0009 

0.0000 

0.0000 

0.0000 

-o.oooo 

8 

n 

4 

0.0006 

-0.0002 

0 . 0000 

-n . nooo 

-o .oooo 

n. onnn 

10 

0 

4 

-0.0001 

0.0002 

0.0000 

-o.oooo 

-n.0000 

— 0 . oooo 

12 

0 

4 

-0.0000 

-0.0001 

0.0000 

0.0000 

0.0000 

-0.0000 
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Table III 

Coefficients in the Series for the Perturbations of Neptune Due to Pluto 


i 

j 

k 

a 

cos 

.10 6 

sin 

4.10 6 

cos 

sin 

cos r - 106 

sin 

0 

u 

0 

0.0361 

0 . 

0 . 

0 . 

- 0.1957 

0 . 

1 

0 

0 

- 2.8317 

- 1.9734 

- 10.7879 

13.4496 

- 0.4121 

0.1763 

2 

0 

0 

7.1479 

0.2280 

0.2795 

- 19.3522 

- 0.5690 

- 0.0708 

3 

0 

0 

- 0.7132 

0.1012 

- 0.0000 

- 0.0000 

- 0.0000 

0.0000 

4 

0 

0 

- 1.8191 

0.9415 

1.6611 

2.9897 

0.2614 

- 0.2765 

3 

0 

0 

- 0.6387 

0.3230 

0.4981 

1.0080 

0.0808 

- 0.0527 

6 

0 

0 

- 0.3069 

0.1337 

0.1731 

0.4327 

0.0608 

- 0.0472 

7 

0 

0 

- 0.1487 

0.0794 

0.0885 

0.1765 

0.0296 

- 0.0324 

,8 

0 

0 

- 0.0709 

0.0467 

0.0580 

0.0799 

0.0146 

- 0.0103 

9 

0 

0 

- 0.0391 

0.0288 

0.0397 

0.0462 

0.0076 

- 0.0066 

10 

0 

0 

- 0.0231 

0.0195 

0.0233 

0.0279 

0.0029 

- O ’. 006 7 

11 

0 

0 

- 0.0132 

0.0128 

0.0132 

0.0155 

0.0021 

- 0.0046 

12 

0 

0 

- 0.0078 

0.0081 

0.0087 

0.0081 

0.0020 

- 0.0027 

13 

0 

0 

- 0.0047 

0.0054 

0.0063 

0.0046 

0.0010 

- 0.0016 

14 

0 

0 

- 0.0027 

0.0037 

0.0043 

0.0029 

0.0003 

- 0.0010 

13 

0 

0 

- 0.0016 

0.0025 

0.0027 

0.0018 

0.0002 

- 0.0008 

16 

0 

0 

-('.0010 

0.0017 

0.0018 

o.oolo 

0.0001 

- 0.0006 

17 

0 

0 

- 0.0006 

0.0012 

0.0013 

0.0005 

0.0001 

- 0 . 0004 

18 

0 

0 

- 0.0003 

0.0008 

0.0009 

0.0003 

0.0000 

- 0.0002 

19 

0 

0 

- 0.0002 

0.0006 

0.0006 

0.0002 

- n.oooo 

- 0.0002 

20 

n 

o 

- 0.0001 

0 . 00^4 

o . nn04 

0.0001 

- o.noon 

- o.oooi 

0 

0 

1 

- 4.0541 

0 . 

0 . 

0 . 

- 0.0027 

o . 

1 

0 

1 

0.0065 

0.0060 

0.0621 

- 0.0524 

0.0003 

- 0.0001 

2 

0 

1 

- 0.0005 

- 0.0002 

0.0034 

0.0046 

0.0004 

0.0002 

3 

0 

1 

1.0487 

3.4693 

6.9393 

- 2.0959 

- 0.3872 

- 0.8137 

4 

0 

1 

0.0024 

- 0.0007 

- 0.0012 

- 0.0042 

- 0.0006 

0.0003 

5 

0 

1 

0.0006 

- 0.0002 

- 0.0004 

- 0.0008 

- 0.0002 

0.0001 

6 

0 

1 

0.0137 

0.0089 

0.0089 

- 0.0138 

- 0.0037 

- 0.0017 

7 

0 

1 

0.0002 

- 0.0001 

- O.COOl 

- 0.0002 

- 0.0000 

0.0000 

8 

0 

1 

0.0001 

- 0.0000 

- 0.0001 

- 0.0001 

- o.cooo 

0.0000 

9 ' 

0 

1 

0.0002 

- 0.0000 

- 0.0000 

- 0.0001 

- 0.0000 

0.0000 

0 

0 

2 

- 0.0007 

0 . 

11.5777 

0 . 

- Q • 0001 

0 • 

1 

0 

2 

- 0.0049 

0.0030 

0.0154 

0 . 0 23 5 

- 0.0001 

- 0.0000 

2 

0 

2 

- 0.0027 

- 0.0016 

- 0.0047 

0.0071 

0.0001 

- 0.0001 

3 

0 

2 

- 0.0003 

- 0.0019 

- 0.0044 

0.0010 

0.0001 

0.0001 

4 

0 

2 

0.0004 

0.0005 

0.0008 

- 0.0006 

- 0.0001 

- 0.0002 

5 

r t 

2 

0.0003 

o.noni 

o. ooo i 

- o .0004 

- o.nooo 

- 0.000 1 

6 

0 

2 

0 . 000 ? 

0.0001 

0.0001 

- 0.0002 

- o • noon 

0.0000 

7 

o 

2 

0.0001 

0.0001 

0.0001 

- 0.0001 

- o.onoo 

0.0000 

8 

0 

2 

0.0001 

0.0000 

0.0001 

- 0.0001 

- 0.0000 

- 0.0000 

0 

0 

3 

0.0010 

0 . 

0.0013 

0 . 

0.0000 

0 . 

1 

n 

3 

- 0.0000 

- 0.0000 

- 0.0001 

0.0000 

- 0.0000 

o. oooo 

3 

o 

3 

0.0005 

- O.0001 

- 0.0003 

- 0.0011 

- 0.0000 

0.0001 

0 

0 

4 

- 0.0001 

0 . 

- 0.0014 

0 . 

0.0000 

0 . 
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Iteration to convergence for the mutual attraction of Pluto and Neptune re- 
quires four iterations and six minutes of machine time on an IBM 7094 "Moon- 
light" system. 

Another case handled by this method of trigonometric expansion is the 
perturbations of the minor planet Hilda due to Jupiter. Again there exists a 
three to two resonance. Table IV gives the elements for the reference orbits 
of the planets. 


Table IV 


Reference Elements of the Planets Hilda and Jupiter 


Planet 

Semi-Major ^ .... 

... Eccentricity 

Axis in a.u. J 

e 

Mean Motion 
in 

Degrees/Day 

Longitude 

of 

Ascending Node 

Argument 

of 

Perihelion 

Inclination 

i 

Mean Anomaly 
at Epoch 
(JD 2433200.5) 


a 

n 

Q 

u 



Hilda 

3.969423 .153760 

?12462789 

228?3400 

49?2610 

7?8460 

245?9780 

Jupiter 

5.203063 .048410 

?08308526 

99?9479 

274?0669 

1?3071 

295?5194 


Table V gives the coefficients of the trigonometric series for a , /3 , and y . 
The unit of time for this series is the period of Hilda divided by 2v . The re- 
duced time unit is necessary because of the large perturbations of Jupiter on 
the minor planet. With a one century time unit, the series coefficients are 
large and convergence is slow. 


Table V 

Coefficients in the Series for the Perturbations of Hilda Due to Jupiter 


i 

j 

k 

a.lU u 

COS sin 

^.10 6 

cos 

sin 

r.io 6 

COS 

sin 

0 

0 

0 

- 120*9943 

0 . 

0 . 

0 . 

- 0.5198 

0 . 

1 

0 

0 

- 306.6840 

- 841.0596 

- 3435.7379 

875.6191 

51.0465 

40. 1660 

2 

0 

0 

1521.6919 

- 2997.8475 

- 8238.6464 - 

4140.0171 

- 46.4376 

1.9010 

3 

0 

0 

- 86.2587 

- 430. 1434 

- 0.0000 

- 0.0000 

o.onoo 

- 0.0000 

4 

0 

0 

727.0608 

718.3318 

1209.8296 - 

1193.0804 

4.4585 

- 62.4435 

5 

0 

0 

- 211.9126 

41 . 605 ? 

4 . 71/6 

346.9106 

4.1683 

19.2708 

6 

0 

0 

228.5693 

- 25.3987 

- 34 . 1563 

- 256.9403 

- 7.1305 

- 7.8846 

7 

0 

0 

- 146.2929 

96.7089 

104.5142 

156.8829 

10.1052 

4.1389 

8 

0 

0 

16.6290 

- 70.0961 

- 75.4031 

- 21.4159 

- 5.1518 

0.7114 

9 

0 

0 

0.0076 

56.6001 

50.9787 

- 3.4319 

2.4647 

- 2.0367 

10 

u 

0 

- 18.2119 

- 33.8711 

- 29.2623 

17.5199 

- 1.1573 

2.2910 

11 

0 

0 

19.7151 

8.1139 

7.0920 

- 18.0331 

- 0.1560 

- 1.5206 

12 

0 

0 

- 15.2318 

- 0.5679 

0.9553 

12.5697 

0.6249 

0.8126 

13 

0 

0 

9.0409 

- 4.8708 

- 4.6395 

- 7.0993 

- 0.6756 

- 0.3321 

14 

0 

0 

- 2.8680 

6.0194 

5.1815 

2.1245 

0.5027 

- 0.0581 

15 

0 

0 

0.0546 

- 4.6413 

- 3.7536 

0.4394 

- 0.2825 

0.2190 

16 

0 

0 

1.6686 

2.7666 

2.1013 

-] .5645 

0.1043 

- 0.2365 

17 

0 

0 

- 2.0506 

- 0.9682 

- 0.6626 

1.7236 

0.0262 

0.1831 

18 

0 

0 

1 .5686 

- 0.0717 

- 0.2072 

- 1.2785 

- 0*0846 

- 0.1046 

19 

0 

0 

- 0.9470 

0.6482 

0.5896 

0.7116 

0.0915 

0.0358 

20 

0 

0 

0.3344 

- 0.7647 

- 0.6330 

- 0 .2205 

- 0.0713 

0.0123 

21 

0 

0 

0.0619 

0.5945 

0.4745 

- 0.0952 

0.0408 

- 0.0347 
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Table V (Continued) 

Coefficients in the Series for the Perturbations of Hilda to Jupiter 


X 

j 

k 

«.10 6 

COS 

sin 

0.1O ( 

cos 

5 

sin 

cos 

sin 

22 

0 

0 

-0.2670 

-0.3515 

-0.2628 

0.2357 

-0.0132 

0.0374 

23 

0 

0 

0.3036 

0.1205 

0.0780 

-0.2481 

-0.0058 

-0.0291 

24 

0 

0 

-0.2362 

0.0360 

0.0435 

0.1864 

C.0148 

0.0166 

25 

0 

0 

0.1379 

-0.1136 

-0.0978 

-0.1026 

' -0.0158 

-0.005 1 

26 

0 

0 

-0.0452 

0.1257 

0.1015 

0.0289 

0.0123 

-0.0027 

27 

0 

0 

-0.0185 

-0.0976 

-0.0761 

0.0199 

-0.0069 

0.0065 

2 8 

0 

0 

0.0494 

0.0562 

0.0416 

-0.0417 

0.0921 

-0.0069 

29 

0 

0 

-0.0536 

-0.0176 

-8.0110 

9.0428 

9.0913 

0.0053 

30 

0 

0 

0.0414 

-0.0C41 

-0.0092 

-9.0320 

-0.0029 

-0.9030 

31 

0 

0 

-0.0236 

C . 0 2 1 8 

0.0162 

0.0173 

0.0030 

0.0008 

32 

0 

0 

Q.0070 

-0.0233 

-9.0185 

“0*0043 

“9.0923 

0.0006 

33 

0 

0 

0.0044 

0.0179 

0.01 37 

-0.0043 

0.0^13 

-0.0013 

34 

0 

0 

-0.00^0 

-0.0101 

-0,0074 

0.0081 

-0.9004 

0.0014 

35 

0 

0 

0.0103 

0.0028 

0.9017 

-8.9981 

“8.8993 

-8.0919 

36 

n 

0 

-9.0079 

9.0021 

0 , 0029 

8 .8969 

8.8896 

9.9906 

37 

0 

0 

0.0044 

-0.0044 

-0.0836 

-9.9032 

-0.8896 

-0.989 1 

38 

0 

0 

-0.0012 

0.0046 

0.0036 

9 .8887 

9.9085 

-9.8801 

39 

0 

0 

-0.0010 

-0.0035 

-0.0O27 

8.8818 

- 0 . 8 9 8 ^ 

8 . 808-3 

40 

o 

0 

9.0029 

0.001.9 

0.8914 

-9,9817 

9.8881 

-9.0093 

41 

0 

0 

-9.0021 

-0.0095 

- o # o o c 3 

9.8816 

9.8.801 

8.9902 

42 

n 

0 

0.0016 

-0.0005 

-0.0005 

-9.0812 

“9.8901 

-0.0001 

4 3 

o 

0 

-9 * 00 n 9 

n . 009Q 

9.0008 

9.8896 

0 • 988 1 

0.9000 

44 

o 

0 

0 • 0092 

-0.0010 

-0.0907 

-0.0901 

“8.0901 

0.0000 

4 5 

0 

o 

0 .0002 

9.0007 

0.0905 

-0.8092 

8.089 ] 

-9.0901 

46 

0 

0 

-0.00O4 

-9.0004 

-0.0903 

9,0094 

-0.8009 

0.0001 

47 

n 

o 

■''.0904 

o . nrn ] 

9.999c 

-8 # 9993 

- 8 . n r 8 o 

-9.0000 

48 

r* 

n 

- r ' . 0 0 9 “3 

O . nno ] 

8 . 9 9 6 2 

0 . 9 9 8 2 

8 # 0888 

8 # 9088 

49 

n 

0 

9.0002 

- 0 . 9 0 o 2 

-9 .0002 

-8.9991 

-0.9909 

-9.0000 

C, 

r. 

n 

-9.009'' 

0 .4092 

9 , PCQ2 

9.9998 

0.9988 

-0 . 0090 

5 1 

n 

0 

-9.00O1 

. 009? 

— 0.9.991 

8 r. 8 8 1 

! * Ul 11 1 

—9 . 8888 

9.0988 

52 

0 

0 

0 . 0 0 o i 

0 • 9 o o l 

0.0801 

-8.9091 

r • n^oo 

-9.0090 

53 


n 

-0.0091 

-4.000C 

-8.0000 

9.00 91 

8.0989 

9.0000 

54 

8 

0 

” .nooi 

# coo c 

— 0 # n on 9 

-8.8891 

-8.8888 

-8 . 98O0 

0 

0 

1 

95.9995 

0. 

0. 

0. 

-4.8773 

0. 

1 

0 

1 

6.7703 

3.8104 

10.9310 

-17.2214 

-0.0497 

0.3242 

2 

0 

1 

-1.0392 

-0.3414 

-0.7731 

7.1091 

-0.1154 

-0.0403 

3 

u 

1 

-564.2345 

9 2 .2646 

16 3. 1807 

1127.3962 

20.0066 

31.7850 

4 

0 

1 

1.1026 

1 .8 002 

-0.2593 

-3.3625 

0.0734 

0.1157 

5 

0 

1 

-2. 1400 

3.3088 

-0.5565 

0.3601 

-0.1196 

0.0147 

6 

J 

1 

1 1 .2623 

-42.9095 

-42.9219 

-11 . j>484 

-2.8424 

0.4044 

7 

0 

1 

-0.6673 

-0.1982 

0.6858 

-0.3163 

-0.0040 

-0.0392 

6 

0 

1 

-U . 7596 

-1.0232 

-0.0838 

0.1660 

0.0110 

-0.0084 

9 

u 

1 

4.0913 

3.1953 

2.1231 

-2.6368 

0.0042 

-0.3227 

10 

o 

1 

0.1613 

-0.1104 

0.0840 

0. 1984 

0.0145 

0.0055 

1 1 

o 

1 

0.32 33 

-0.1223 

-0.0931 

-8.0432 

-0.0045 

0.0037 

12 

0 

1 

-0.6166 

0.3736 

0.2 38 7 

0.3072 

O.0385 

0.0254 

13 

o 

1 

- 0 . 0 1 3 3 

9 * 0 3 9 6 

-8.0604 

9.9286 

-9.9928 

9. 0049 

14 

0 

1 

0.0167 

0.0825 

0.0194 

-0.0393 

- 0 • 0 0 1 0 

-0.0032 

15 

r: 

1 

-0.0332 

-0.1003 

-0.0402 

9.0349 

-0.0045 

0.0050 

16 

0 

1 

-0.0010 

-O.Q186 

-0.0123 

-0.0219 

“9*0019 

-0.0011 

17 

o 

1 

-0.0142 

9! .0067 

0.0166 

0.0080 

0.0016 

-0.0003 

18 


1 

0 . 0 1 3 1 

-0.0079 

-0.0118 

-C.OOhI 

-0.0010 

-0.0003 

19 

'■) 

1 

0.0102 

0.0051 

0.0088 

-0.0058 

0.0004 

-0.0009 

20 

0 

1 

-9.00 80 

-0.0047 

-0.0033 

0.0074 

0.9881 

0.0808 

21 


1 

0.0040 

O.00O7 

-o.ono? 

-8.9955 

-8.0992 

-0.0004 

22 

0 

1 

-9.0037 

0.0045 

0.0028 

9.0038 

0.9884 

0.0002 

2 3 

i; 

1 

9.0015 

-0.0032 

-0.0034 

-0.0014 

-o . 9994 

0.0001 

24 

0 

1 

9.0094 

0.0029 

0.9927 

-9 ,9004 

8.0992 

“0.0002 

2 6 

0 

1 

-0.0020 

-0.0029 

-0.0017 

0.0014 

-0.0091 

0.0002 

26 

v 

1 

0.0018 

0.0006 

0.0806 

-0.0016 

-0.0000 

-0.0002 

27 


1 

- >.0015 

9.0003 

9.9002 

8.9913 

8.9081 

o.nooi 

28 


1 

9.0019 

-C .0009 

-0.0907 

-9.0988 

-0.9901 

-0.0000 
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Table V (Continued) 

Coefficients in the Series for the Perturbations of Hilda Due to Jupiter 


«.io 6 4-io 6 y.io 6 


i 

j 

k 

COS 

sin 

cos 

sin 

cos 

sin 

29 

/-> 

1 

-P.0003 

0.0099 

0 . 09P8 

9.9007 

o.oooi 

-9.0009 

30 

0 

1 

-0.0002 

-0.0008 

-0.0006 

0.0001 

-0.0091 

0.0000 

31 

0 

1 

0.0004 

0.0005 

0.0004 

-0.0003 

0.0000 

-9.0001 

32 

0 

1 

-0.0005 

-9.0001 

-0.0001 

9 .9094 

0 . 0009 

0 . 0090 

33 

q 

1 

9.0094 

-9.0091 

- 0.0001 

-p . POO 3 

- 9 . 0 0 O P 

-O.0000 

34 

0 

1 

-0.0002 

0.0002 

0.0002 

9.0002 

0.9000 

0.0000 

3 5 

0 

1 

0.0001 

-0.0007 

-0.0002 

-o.ooon 

-o.oono 

0.0000 

36 

0 

1 

0 . 0000 

0.0997 

p . pop l 

-p . 0000 

9.0009 

-9.0900 

37 

o 

1 

-o . oon i 

-9.9001 

-o.OOOl 

9 . n 0 9 1 

-P .9090 

9.0009 

38 

0 

1 

o.OOOl 

0.0009 

0.0000 

-9.009 1 

- 0.0000 

-0.0900 

39 

o 

1 

-0.0001 

o.ooon 

0. 009P 

9.0091 

0 . 0 0 or 

0.0000 

^0 

n 

1 

n .OOP 1 

-p.nooi 

-o . Ppo ^ 

-n . pppo 

— p . n r « 0 

- P , POP p 

41 

0 

1 

-0.0000 

9 . 0 09 1 

0.0000 

9.9099 

P . OPOP 

-p.ooon 

0 

0 

2 

-0.5686 

0. 

-53.7978 

0. 

-0.0252 

0 . 

1 

0 

2 

-0.0831 

0.0701 

0.6767 

-0.5963 

9.0101 

-0.0142 

2 

0 

2 

-0.1962 

-0.1008 

0.4002 

-0.4359 

0.0017 

0.0143 

3 

0 

2 

-0.1275 

0.1611 

0.6963 

0.4037 

-0.0232 

0.0017 

4 

0 

2 

-0. 0434 

-9.0090 

0. 1174 

9.1888 

0 . 0922 

0.0007 

5 

0 

2 

0.1065 

-0.0641 

-0.0652 

-0.0 1 36 

-0.0018 

0.0006 

6 

0 

2 

0. 1733 

-0.0209 

0.0438 

9.9957 

0.9078 

9.0058 

7 

0 

2 

-0.0016 

0.0016 

-0.0329 

0.0159 

-0.0003 

0.0011 

8 

0 

2 

0.0124 

0.0327 

0.0064 

-0.0181 

-0.0003 

-0.0008 

9 

0 

2 

-0.0219 

0.0678 

-0.0129 

0.0090 

-0.0018 

9.0012 

10 

0 

2 

-9.0004 

-9.0048 

-0.004 1 

-9.n 0 87 

-0.0993 

-9.00 02 

11 

0 

2 

-0.0080 

0.0019 

0.0059 

9.0025 

0 . 0004 

“9.900 1 

12 

u 

2 

-0.0164 

-0.0148 

-0.0032 

-0.0015 

-0.0002 

-0.0003 

13 

0 

2 

0.0029 

0.0019 

0.0028 

-O.0015 

0.000] 

-0.0002 

14 

0 

2 

-0.0009 

-0.0018 

-0. onl 0 

0.9022 

p . nnoo 

0 . 0002 

15 

0 

2 

9.0059 

-P.0028 

0 . 0991 

-9.0914 

-9.9099 

-O.OOOl 

16 

0 

2 

-0.0009 

0.001 3 

0.000 7 

9.9011 

9.0901 

0.9000 

1^ 

n 

2 

0.00P5 

-0.0096 

-0. 0900 

-p . 0 004 

-0.099] 

0 . 0 0 p 0 

18 

0 

2 

0.0003 

0.9070 

0.0007 

-0.0009 

0 • 0 9 0 0 

- O . 9 0 0 0 

19 

0 

2 

-0.0095 

-0.0005 

-0.0004 

0 . 0 OP 3 

- 0 . 9909 

9.0000 

29 

p 

2 

''.0904 

9.9007 

0 . nop? 

-p . n 004 

— 0 . O 0 P P 

- 0 . OOP p 

21 

0 

2 

-0.0006 

-o.ooon 

0 . noon 

9 . 9 903 

P . OOPP 

p . OPPO 

22 

0 

2 

0.0093 

- 0.0002 

- 0.0001 

• -P .9992 

-0 .OOOP 

-0.0099 

29 

p 

2 

-'’.onoi 

p . 0 n 0 7 

n . ~pp? 

9 . 9 99 1 

p . A p p p 

- P . PPPP 

24 

0 

2 

0.0000 

- 0.0002 

- 0.0001 

0.0099 

- 0.0000 

o.ooon 

26 

0 

2 

0.0001 

0.0001 

0.0001 

-9.0001 

0.0000 

- 0.0000 

26 

0 

2 

- 0.0001 

- 9 .noon 

- 0 , npnn 

9 . p 90 1 

0 . 0 9 0 0 

9. 0000 

27 

0 

2 

0.0091 

0.0000 

-0 . 0009 

-9 .9091 

-0 • oopp 

-9 . 0009 

28 

0 

2 

-0 .000-1 

o.ooon 

0.009 0 

9 .0999 

0 . 09PP 

9.0009 

0 

u 

3 

0.0014 

0 . 

0.2027 

0. 

0. 0002 

0 . 

1 

n 

3 

0.0003 

-9.0006 

0.0046 

9.0076 

-r . noon 

0.0001 

2 


3 

9.0006 

9 .0007 

-n. '■'026 

9.0096 

q . pppp 

— 9 . 0 0 P 1 

3 

o 

3 

0.0045 

9.0201 

-0.0597 

0.0003 

-o.ori 3 

9.0006 

4 

0 

3 

—9.0091 

p .poop 

-0 . POPP 

-0 . 0 opq 

_p . pppp 

-p. pppp 

5 

n 

3 

-9.0008 

-9.009] 

0 . 0 0 p 3 

— 9 . 9 9 P 2 

r .poop 

-p . PPPO 

6 

0 

3 

-0.0109 

-0.0037 

0 . CO 1 2 

-0.0033 

- 0.0000 

-0.0002 

7 

0 

3 

O'. 000 1 

-0.0002 

0.0002 

o.OOOl 

0.0 OOP 

-p.onoo 

8 

0 

3 

0.009-7 

-O.P092 

- 0 , nnoo 

9 . OOP 1 

0 . POPP 

P . 0000 

9 

0 

3 

9.007 3 

-0.9074 

0.0999 

p .9993 

P . OPPO 

p. O opp 

10 

p 

3 

o.ooon 

P.OPni 

-0.0001 

P.090 1 

P . POPO 

0 . noon 

1 1 

n 

3 

9,9090 

o .non } 

-0 . oppp 

-9 . P PPP 

-P . OPPP 

-P. OOOP 

1? 

q 

3 

0.0003 

9.0098 

-0 . OOpp 

P . P PPP 

— p . p p P P 

P . POPP 

15 

0 

3 
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